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Abstract 

Doubly special relativity (DSR) is an effective model for encoding quantum gravity in flat space- 
time. As a result of the nonlinearity of the Lorentz transformation, the energy-momentum dis¬ 
persion relation is modified. One simple way to import DSR to curved spacetime is “Gravity’s 
rainbow”, where the spacetime background felt by a test particle would depend on its energy. Fo¬ 
cusing on the “Amelino-Camelia dispersion relation” which is = mf — g{E/mp)^] with 

n > 0, we investigate the thermodynamical properties of a Schwarzschild black hole and a static 
uncharged black string for all possible values of g and n in the framework of rainbow gravity. It 
shows that there are non-vanishing minimum masses for these two black holes in the cases with 
g < 0 and n > 2. Considering effects of rainbow gravity on both the Hawking temperature and 
radius of the event horizon, we use the geometric optics approximation to compute luminosities 
of a 2D black hole, a Schwarzschild one and a static uncharged black string. It is found that the 
luminosities can be significantly suppressed or boosted depending on the values of g and n. 
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I. INTRODUCTION 


The Hawking radiation was a remarkable prediction of quantum held theory in curved 
spacetime. As Stephen Hawking demonstrated[l|, a Schwarzschild black hole emits radia¬ 
tion just like an ordinary blackbody at temperature T = hK/2kB7i, where n is the surface 
gravity. Soon after this discovery, it was realized that there might be the ’’trans-Planckian” 
problem 2j . It appears that the Hawking radiation originates from the modes with huge 
initial frequencies, well beyond the Planck mass nip, which undergo exponential high grav¬ 
itational red-shifting near the horizon. As a result, the Hawking’s prediction relies on the 
validity of quantum held theory in curved spacetime to arbitrary high energies. On the other 
hand, quantum held theory is considered more like an ehective held theory of an underlying 
theory whose nature remains unknown. This observation poses the question of whether any 
unknown physics at the Planck scale could strongly inhuence the Hawking radiation. 
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Although a complete understanding of the trans-Planckian problem requires a full the¬ 
ory of quantum gravity, there are various attempts using effective models to address this 
problem. Among them is Doubly Special Relativity (DSR)j3-6|, where the transformation 
laws of special relativity are modified at very high energies. In DSR, the energy-momentum 
dispersion relation for a massive particle of mass m is modified to 


E^r{E/m,)-p^g^ {Elm,)=m\ 


( 1 ) 


where rup is the Planck mass, and / {x) and g (x) are two unknown functions with the 
following properties 

lim f (x) = 1 and lim g (x) = 1. (2) 

x^O x^O 

It has been shown that the modified dispersion relation (MDR) might play a role in astro¬ 
nomical and cosmological observations, such as the threshold anomalies of ultra high energy 


cosmic rays and TeV photonsj7-ll2|. Moreover, thermodynamics of black holes have been 


explored in the framework of the MDR 
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18|. 


One of the most popular choice for the functions / (x) and g (x) has been proposed by 


Amelino-Camelia et al. 


19 


20l |. which gives 


f(x) = l and g (x) = y/T^ 




( 3 ) 


Usually one has n > 0. As shown in 2^, this formula is compatible with some of the 
results obtained in the Loop-Quantum-Gravity approach and reflects the results obtained in 
K-Minkowski and other noncommutative spacetimes. Phenomenological implications of this 
“Amelino-Camelia (AC) dispersion relation” are also reviewed in 20|. 

Nevertheless, the non-linear realization of the Lorentz transformation in the framework of 
DSR results in a very complicated definition of the dual position space. To circumvent this 


difficulty, Magueijo and Smolin 


2l| proposed the “Gravity’s rainbow”, where the spacetime 


background felt by a test particle would depend on its energy. Consequently, the energy 
of the test particle deforms the background geometry and hence the dispersion relation. 
As regards the metric, it would be replaced by a one parameter family of metrics which 
depends on the energy of the test particle, forming a “rainbow metric”. Specifically, for the 
Schwarzschild solution, the corresponding “rainbow metric” is 

{ JT' /im T 

^pj 


ds^ = 1 


2GM\ dt"^ 1 

r ) P {E/nip) g^ {E/rrip) 


2GM 


+ P + sin^ Odp') 


( 4 ) 
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which is the spherically symmetric solution to the distorted Einstein’s Field equations 


given m 


2l|. Note that E is the energy of the probing particle. To obtain the modi¬ 
fied Schwarzschild metric from the usual one, it appears from eqn. (jTj) that we can simply 
make replacements dt —)■ dt/ f [E/rrip) for the time coordinate and dx* —)■ dx'^/g {E/rrip) for 
all spatial coordinates. In fact, such procedure also works for the black objects besides the 
Schwarzschild metric j^. The rainbow gravity formalism has received a lot of attentions 


recently, for instance, in cosmology 


23 


26 1 and black hole physics 27-34| 


In this paper, we will study thermodynamics and luminosities of black holes in the frame¬ 
work of rainbow gravity. The remainder of our paper is organized as follows. In section [ITl 
the deformed Hamilton-Jacobi equations for scalars, spin 1/2 fermions and vector bosons are 
derived in the framework of rainbow gravity. We then solve the Hamilton-Jacobi equations 
to obtain tunneling rates. The temperatures and entropies of a rainbow Schwarzschild black 
hole and a rainbow static uncharged black string is computed in section HTTl In section HVl we 
calculate the luminosities of a 2D rainbow black hole, a 4D rainbow spherically symmetric 
one, and a 4D rainbow cylindrically symmetric one. Section |V] is devoted to our conclusion. 
Throughout the paper we take Geometrized units c = G = 1, where the Planck constant h 
is square of the Planck mass m^. 

II. HAMILTON-JACOBI METHOD 


After the Hawking’s original derivation, there have been some other methods proposed to 
understand the Hawking radiation. Recently, a semiclassical method of modeling Hawking 


radiation as a tunneling process has been developec 


method was first proposed by Kraus and Wilczek 


35 


and attracted a lot of attention. This 


36[ |. which is known as the null geodesic 


method. They employed the dynamical geometry approach to calculate the imaginary part 
of the action for the tunneling process of s-wave emission across the horizon and related it 
to the Hawking temperature. Later, the tunneling behaviors of particles were investigated 
using the Hamilton-Jacobi method 


37 


39l |. In the Hamilton-Jacobi method, one ignores 


the self-gravitation of emitted particles and assumes that its action satisfies the relativistic 
Hamilton-Jacobi equation. The tunneling probability for the classically forbidden trajectory 
from inside to outside the horizon is obtained by using the Hamilton-Jacobi equation to 
calculate the imaginary part of the action for the tunneling process. 
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In this section, the Hamilton-Jacobi equations for scalars, spin 1/2 fermions and vector 
bosons in the rainbow metric are derived. The particles’ tunneling rates across the event 
horizon r = oi the rainbow metric ([S]) are then computed by solving the Hamilton-Jacobi 
equations. We here consider a static black hole with the line element 


ds = B (r) dt — 


B(r] 


— C (r^) hab {x) dx'^dx^, 


( 5 ) 


whose rainbow metric is obtained by dt —)• dt/f {E/irip) ,dr dr/g{E/mp), and dxi —>■ 
dxi/g {E/mp) 


ds^ = 


B (r) dt^ 


dr"^ 


C (r^) hab (x) dx°‘dx'^ 


( 6 ) 


P {E/mp) g^ {E/mp) B (r) g^ {E/mp) 

The function B (r) has a simple zero at r = r/j with B' {r^) being hnite and nonzero. The 
vanishing of B (r) at point r = Vh indicates the presence of an event horizon. For simplicity, 
we assume that the particles are massless. 

In the rainbow metric ds^ = g^i, {E) dx^dx'^, the massless scalar held (j) obeys the Klein- 
Gordon equation 

WV,/) = 0, (7) 

where is the covariant derivative associated with g^y {E) and the index /i is lowered or 
raised by g^i^ {E). Making the ansatz for 0 which is 


= exp 


a 

'hi ’ 


( 8 ) 


substituting it into eqn. o. and expanding eqn. (JTj) in powers of h, the leading order gives 
the Hamilton-Jacobi equation for a massless scalar particle 


~g^, {E) aWJ = 0. 

The Dirac equation for a spin-1/2 fermion held 'i/ takes the form of 

V- = 0, 


( 9 ) 


( 10 ) 


where = ^u/^’^Eab, and U!//° are the Lorentz spinor generator and spin connection in 
the rainbow metric, respectively, and = ‘^gfj.u{E). To obtain the Hamilton-Jacobi 

equation for the fermion, the ansatz for is assumed as 

il 


= exp ( -^ ) u. 
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( 11 ) 






where f is a slowly varying spinor amplitude. Substituting eqn. (ITT]) into eqn. OB, we find 
to the lowest order of h 

= 0 . ( 12 ) 

Multiplying both sides of eqn. (IT^ from the left by one gets 

[g^AE)d^Id‘'I]v = 0. (13) 

Since v is nonzero, the Hamilton-Jacobi equation for a massless fermionic particle is also 
given by eqn. (E]). 

The Maxwell’s equations for a massless vector field is 

= 0, (14) 

where = V We then make the WKB ansatz 

= O/, exp , (15) 

where is the polarization vector and I is the action. Plugging the WKB ansatz into eqn. 
oa, we find that leading order of h gives 


{E) {a.djd^l - a^djdj) = 0. 


(16) 


To simplify eqn. flTHl) . one could impose the Lorentz gauge in the curved spacetime 


= 0 , 


(17) 


where the leading order is 


g-’{E) aAI = 0. 


(18) 


By plugging eqn. f[T5]l into eqn. fITB]) . it shows that the Hamilton-Jacobi equation for a 
massless vector boson is also given by eqn. ([9]). 

From eqn. (j9]) , one finds that the Hamilton-Jacobi equation for a massless particle in the 
rainbow metric ([6]) becomes 


P {E/nip 


jdjy 

B(r) 


= / {E/mp) 


Df ^ 70 n2 , {x)daldbl 
B{r){drl) +-- 


(19) 


To solve the Hamilton-Jacobi equation for the action J, we can employ the following ansatz 


/ = -Et + JF (r) + 0 (x), 


( 20 ) 
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where E is the particle’s energy. Plugging the ansatz into eqn. (IT^ . we have differential 
equations for W (r) and 0 {x) 


(x) daQ {x) dbQ (x) = A, 
drW± (r) =Pr = 


_^±_ ±1 lE^PiE/mp) 


( 21 ) 


B (r) y {E/mp) C {r^) ’ 

where +/— denotes the outgoing/ingoing solutions and A is a constant. Using the residue 
theory for the semi circle, we get 

±TTEf {E/rrip) 


Im W± (r) = 


2k g{E/mp) 


( 22 ) 


where k = As shown in 


40|, the probability of a particle tunneling from inside to 
2 


h 


(Im 1U+- Im IU_) 


outside the horizon is 

Pemit OC exp 

There is a Boltzmann factor in Pemit with an effective temperature, which is 

g{E/mp) 


Teff - To 


where we dehne To = fy and take ks = 1- 




(23) 


(24) 


III. THERMODYNAMICS OF BLACK HOLES IN RAINBOW GRAVITY 

We now use the Heisenberg uncertainty principle to estimate the black hole’s temperature. 
For a massless particle, the modihed dispersion relation ([1]) becomes 

xf (x) 




= y, 


(25) 


where x = E/irtp > 0 and y = p/mp > 0. To obtain the black hole’s temperature from eqn. 
fl2Tll , eqn. (|25ll is needed to be solved for x in terms of y. In fact, one could have for x 


X = yh (y), 


(26) 


is inverted to get the function h {y) and lim h{y) = 1. The Heisenberg 


where eqn. 

uncertainty principle gives a relation between the momentum p of an emitted particle and 
the event horizon radius of the black hole 1^| 


y = p/rrip ~ Sp/rUp ~ h/mp5x ~ mp/vh- 


(27) 
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Substituting eqn. IHT} into eqn. (1^ . we have for the black hole’s temperature 

T = To- = Toh('^y (28) 

P \rh J 

The range of the left-hand side(LHS) of eqn. fl2^ determines the ranges of the values of Vh- 
Specihcally, the maximum value of the LHS of eqn. fl25|) . which is denoted by j/cr, gives that 
If Vcr is hnite, it always predicts the existence of the black hole’s remnant. For some 
functions / (x) and g (x), the domain of the LHS of eqn. (|25ll might be [0, Xcr]/[0, Xcr) with 
Xcr being hnite rather than [0, cxo). Thus, it gives that the energy of the particle E < rUpXcr- 
If the domain is [ 0 , cxo), we simply set Xcr = oo. 

For the AC dispersion relation given in eqn. ([3]), eqn. fl 2 ^ becomes 


X 

"v/r^-^Tx” 


= y- 


(29) 


If r; > 0, we hnd that ycr = 0. However, there is an upper bound Xcr = on x to make 

the LHS of eqn. real. If r; < 0, we have Xcr = oo and ?/cr = cxo for 0 < n < 2 and 

Xcr = oo and i/cr = (—lor n = 2. For the case with rj < 0 and n > 2, the LHS of eqn. 
fl29|) has a global maximum value yo at Xq, where we dehne 


Xo = 


yo = 



(30) 


Thus, it would appear that y < yo and x < oo since x can go to inhnity. In FIG. [T], we plot 
an example with g = —1 and n = A. If one solves eqn. (l2^ for y < yo in terms of x, there 
are always two solutions xs and xl, where xs < xo < xl- However, only the solution xs is 
Taylor expandable in r]. The solution x^ is a ’’runaways” solution since it does not exist in 
the limit of 77 —)■ 0. In it has been argued that ’’runaways” solutions are not physical 
and hence should be discarded. A similar argument was also given in the framework of the 


Generalized Uncertainty Principle in 


44l |. Therefore, we will discard ’’runaways” solutions 


and keep only the solution xs in this paper. In this case, we have Xcr = xo instead of 
Xcr = 00 . We list Xcr and ycr for various choices of n and rj in TABLE [B 

To calculate the black hole’s temperature from eqn. fl28l) , we need to solve eqn. fl2^ for 
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h {y). Nevertheless, it is non-trivial to get h {y) for a general n. For n = 1, 2,4, we have 


h{y) = < 


-rjy + \/4: + rfy"^ 


(1 + yy 


2N-1/2 


y 


-2 


\/l+4?jy4_i 

2r, 


If 1 / <C 1, one has x <C 1 and hence eqn. (12^ becomes 


n = 1 
n = 2 

n = 3 




(31) 


(32) 


which gives 

h{y) = 1-^ + 0 (y^") . (33) 

Since x = E/nip and y = p/nip, the AC dispersion relation for a massless particle fl2^ 
leads to the group velocity of the particle 


Vg dx a/I - 2 (1 - 

If ?7 > 0, eqn. fIMl) gives Xg < 1, which means the particle is subluminal. If rj < 0, Xg > 1 
and hence the particle is superluminal. 

We now calculate the temperatures and entropies of a Schwarzschild black hole and a 
static uncharged black string with negative cosmological constant. 



Xcr 

Vcr 


/nip 


/nip 

Lines in figures 

rj = 0 

oo 

OO 

0 

oo 

0 

0 

Blue Solid 

y > 0 

Pj-l/n 

oo 

0 

g-l/n 

An 

0 

0 

Black Solid 

y < 0,0 < n < 1 

OO 

oo 

0 

oo 

0 

oo 


y < 0,n = 1 

oo 

oo 

0 

oo 

0 

Sa^nip 

y 4n 

Black Dashed 

y < 0,1 < n < 2 

oo 

oo 

0 

oo 

0 

0 


y < 0,n = 2 

oo 


nip 

oo 

a^rrip 

oo 

Red Dashed 

y < 0,n > 2 

Xo 

yo 

nip 

‘^yo 

XQ 

An 

a^rrip 

^y'o 

^ 

Red Solid 


TABLE I: The values of Xcr, Vcr, /nip, and /nip for various values of {n,y). 

The superscripts SC and BS stand for the Schwarzschild black hole and black string, respectively. 
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FIG. 1: Plot of the curve of^=== with r? = — 1 and n = 4. The curve has a maximum value 
2/0 = ^at xo = 1. For y < yo, there are xs and xl with xs < xl satisfying = y- The 

solution xi, is considered runaways since it does not exist in the limit of ?/ —>■ 0. Thus, only the 
solution xs is considered in the paper. 

A. Schwarzschild Black Hole 


For a Schwarzschild black hole with mass M, one has B (r) = 1 - — and rh = 2M. 

Thns, the temperatnre is 

= (35) 

SttM \2MJ ^ ’ 

If M ;:§> irip, we have from eqn. (l33il that 


jiSC ^ '^p 


StiM 


1 - 


T] m. 
2»^+i M"- 


p +0/ < 




(36) 


The minimnm mass of the black hole is given by 


MSC = 


2|/c 


(37) 


When the mass M reaches , the temperatnre of the black hole is denoted by . Eqn. 
gives that 

(38) 


rjiSC _ 


For r] < 0 and n > 2, //cr is hnite and hence the black hole shonld have non-vanishing 
minimum mass . This implies an existence of black hole remnant due to rainbow 
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FIG. 2: Plot of the temperature /rup against the mass M/mp for a Schwarzschild black hole. 
All the lines asymptotically approach = 0 as M/mp —>■ oo. The blue line is usual case, where 
blows up as M ^ 0. The red dot is the end of the red solid line, where the black hole has a 
remnant = l/\/2mp. In this case, does not blow up as M —. The black dotted 

line is the asymptotic line of the red dashed line as M ^ = 0.5mp, which is the black hole’s 

remnant. In this case, blows up as M —)• . 


gravity. By eqn. (I38|) . we find that Tj!*" for n = 2 is infinite while Tj!*" for n > 2 is 
which is hnite. For r/ < 0 and 0 < n < 2, we hnd that = 0 and = oo. In this 
case, the black hole would evaporate completely while its temperature increases and hnally 
becomes inhnity during evaporation, just like the standard Hawking radiation. For rj > 0, 
the black hole would also evaporate completely. However, the temperature of the black hole 


IS 


a hnite value - 


-l/rij 


An 


at the end of the evaporation process. We list and for 

all the possible values of r] and n in TABLE [H In FIG. [21 we plot the temperature /nip 

against the black hole mass M/mp(both in Planck units), for examples with {T],n) = (1,2), 
{ri,n) = ( — 1,1), {r],n) = (—1,2), and {ri,n) = (—1,4). The standard Hawking radiation is 
also plotted as a blue line in FIG. El 

From the hrst law of black hole thermodynamics dS = dM/T^^, the entropy associated 
with the temperature fl35ll can be given by 






<msc 


dM 

J'SC 


^Vcr 


= 2ti 


dy 


y'^h{yy 


(39) 
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FIG. 3: Plot of the entropy S against the 
mass M/rUp for a Schwarzschild black hole. 


FIG. 4: Plot of the entropy of a Plank length 
snip against the mass per unit length M for 
a static uncharged black string. 


where we use ycr = rj = 0, we have h{y) = 1 and ycr = oo. Thus, eqn. 

the Bekenstein-Hawking entropy 

AttM^ _ A 
rrii Ah 


S = 


gives 


(40) 


where A = is the area of the black hole. If M ;§> Trip {A^ h)^ eqn. 039 p gives the 
entropy up to the subleading term 




2 — n 

+ S (ifl) ^ 0 < n < 2 




A + M In ^ 

2 iwh 


ih 


+ Sn 


n = 2 
n > 2 


(41) 


where Sq is a constant and we use eqn. fl33p . The leading terms of eqn. flTip are the familiar 
Bekenstein-Hawking entropy. For n = 2, we obtain the logarithmic subleading term. In 
FIG. 131 we plot the entropy S against the black hole mass M/rUp, for examples with r; = 0, 
(h, n) =^( 1 , 2 ), (p, n) = (- 1 , 1 ), ( 77 , n) = (- 1 , 2 ), and ( 77 , n) = (- 1 ,4). 

In [3^ , the temperature and entropy of a rainbow Schwarzschild black hole has also been 
computed in the case with 17 > 0 and n = 2. The temperature and logarithmic subleading 
term of the entropy obtained in our paper are the same as those i n |34l| . By contrast, the 
temperature and entropy were calculated in the case with 77 > 0 in 33|, where it was found 


that the black hole had a non-vanishin 
the minimum mass. As pointed out in 


; minimum mass and the temperature was zero at 


and 


34| . the difference between the results in our paper 


34| and those in 33(] stems from the fact the MDR was partially used in SSj. In fact. 
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the author in 331 used the MDR to get Teff while the ordinary dispersion relation rather 
than the MDR was used to obtain the relation between the energy E of an emitted particle 
and the event horizon radius r^. 


B. Static Uncharged Black String 


We consider the Einstein-Hilbert action in four dimensions with a negative cosmological 
constant A, which is given by 

1 




IbTT 


j d‘x^f^(R-2h), 


(42) 


where g is the determinant of the metric and R the Ricci scalar. If we assume that the 
solution of the Einstein equations is cylindrically symmetric and time-independent, the line 
element for a static uncharged black string becomes 4^ 

h \ , o / o o b 


ds = a r 


, 2 2 


where 


dt^ — ( — 

ar J V oir 


A 


-1 


dr'^ — [dO'^ + a^dz^') , 


(43) 


a‘ = h = 4M 
o 

-cx) <f<oo, 0<r<oo, —cxo<z<cx), 0<6 *< 27r. 


(44) 


Here, M is the mass per unit length of the 2 ; line of black string. 

For the static uncharged black string, we have B (r) = ~ /ot^ and 

K = Thus, eqn. fl28|l leads to the black string temperature 


Sarripb^^^. /amp 
- = - - - h ' 


dvr V 6^/3 

If 6^/^ S> amp (M 3> a^m^j , we have for the temperature 

7] a^mp 


dvr 


1 - 


2^+1 5 n /3 


+ 0 




l)2n/3 


(45) 


(46) 


If n > 1, one has 


-iBS 


—)■ To as M —)■ 00 , 


where Tn = 




47r 


is the black string temperature in the usual case. Since r/j > ^, the 


minimum mass per unit length of the black string is given by 


„„ a^ml 
MBS = - P 

cr 4„,3 ■ 




(47) 
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FIG. 5: Plot of the temperature /nip against the mass per unit length M for a static uncharged 
black string. The lines with n > 1 asymptotically approach the blue line as M —?■ oo. The blue line 
is the usual case, where —>• 0 as M —)• 0. The red dot is the end of the red solid line, where 
the black string has a remnant = l/\/2. The black dotted line is the asymptotic line of the 

red dashed line as M ^ = 0.25, which is the black string’s remnant. In this case, blows 

up as M —>• . 


When M reaches , the temperature of the black string becomes 


Tc/ _ x^r 


nip dvr y, 


(48) 


For the usual case with 77 = 0, we have h{y) = 1 and ycr = 00 . The minimum mass per 
unit length = 0 and goes to zero as M —)■ 0. For 77 < 0 and n > 2, y^r is hnite 
and hence the black string should have non-vanishing minimum mass per unit length , 
which indicates an existence of black string remnant. Using eqn. (1481) . we can show that 
for 77 = 2 is inhnity while for 77 > 2 is hnite. For 77 > 0, the black string would evaporate 
completely with goes to zero as M —)• 0, which is the same as in the usual case. For 
?7 < 0 and 0 < tt. < 2, the minimum mass per unit length is zero. Nevertheless, we 

hnd that = 00 for 0 < 77 < 1, for 77 = 1, and = 0 for 0. We list 

and for various values of 77 and n in TABLE [H In FIG. [5] where we set anip = 1, 
we plot the temperature /nip against the mass per unit length M, for examples with 
77 = 0, ( 77 , 77 ) = (1, 2), ( 77 , 77 ) = (-1,1), ( 77 , 77 ) = (-1, 2), and ( 77 , 77 ) = (-1,4). 

Using the hrst law of black hole thermodynamics, we have for the entropy per unit length 
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associated with the temperature fHSD 


s = 


rM 


I msc 


dM Trarrip 


dy 


J'SC 


rrir 


I y^h {y)' 

If M S> 1, eqn. gives the entropy per unit length up to the sub leading term 


(49) 


na I anri fe(2 "-(/S n ^ ^ ^ O 

"o--T “r TTTTi -7 7-777777 U \ / 6 \ Zt 

2 2{2-n) {anipY 


S ~ < 


_j_ narj 

'»/2,v-vi2 O ryu/tyi 


2 


amp 


TTQ _|_ 

2 a2^2 *0 


n = 2 
n > 2 


(50) 


where sq is a constant. In FIG. [3l we plot rripS against M, for examples with rj = 0, 
iVP) = (1,2), {T],n) = (-1,1), {v,n) = (-1,2), and {r],n) = (-1,4). 

IV. LUMINOSITIES OF BLACK HOLES IN RAINBOW GRAVITY 

For particles emitted in a wave mode labelled by energy E and quantum numbers i, we 
hnd that 

(Probability for a black hole to emit a particle in this mode) 

E \ 

X (Probability for a black hole to absorb a particle in the same mode). 


= exp 


Z 


eff 


where Tg// is given by eqn. 


obtained by Hartle and Hawking 


The above relation for the usual dispersion relation was 


46| using semiclassical analysis. Neglecting back-reaction. 


exp , „ 

^ ' Eff 

mode with E and i 


where we dehne 


detailed balance condition requires that the ratio of the probability of having N particles 
in a particular mode to the probability of having N — 1 particles in the same mode is 
^ ' . One then follows the argument in ^ to get the average number nE,i in the 

nE,i = n , ( 51 ) 

n{x) = -(52) 


exp X — (—1)^ 

Note that e = 0 for bosons and e = 1 for fermions. 

As discussed in section mu there is an upper bound rUpXcr on the energy E of the particle. 
Moreover, another upper bound comes from the requirement that nothing can be emitted 
that lowers the energy below the remnant mass of a black hole. Thus, one has E < M — Mcr, 
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where M is the mass of the black hole and Mcr is the remnant mass. Considering both npper 
bonnds, we have for E 


E "E: f^max — min ^TTlpXQ^f^ . 


(53) 


In 


40l . 1471148| . we have fonnd that the total Inminosity in the MDR case is given by 


L 




(54) 


where E is the energy of the particle, i are qnantnm nnmbers needed to specify a mode 
besides E, and |Tj {E)f is the greybody factor. Usnally, Tj (E) represents the transmission 
coefficient of the black hole barrier which in general can depend on E and i. The relevant 
radiations usnally have the energy of order hM~^ for a black hole with the mass M and 
one hence needs to use the wave equations to compute Ti (E). However, solving the wave 
equations for Tj (E) could be complex. On the other hand, we can use the geometric optics 
approximation to estimate \Ti{E)f. In the geometric optics approximation, we assume 
E ^ M and high energy waves will be absorbed unless they are aimed away from the 
black hole. Hence |Tj {E)^ = 1 for all the classically allowed E and i, while |Tj {E)\^ = 0 
otherwise. For the usual dispersion relation, the well-known Stefan’s law for black holes is 
obtained in this approximation. 

In the remaining of the section, we will use the geometric optics approximation to calcu¬ 
late luminosities of a 2D rainbow black hole, a 4D rainbow spherically symmetric one, and 
a 4D rainbow cylindrically symmetric one. For simplicity, we assume that the particles are 
massless. 


A. 2D Black Hole 


Suppose the metric of a 2D black hole is given by 


ds^ = B (r) dt^ - —dr'^, 

B (r) 


(55) 


where B (r) has a simple zero at r = Vh- In this case, we have |T {E)f = 1 when E < E„ 
and |T {E)f = 0 otherwise. By eqn. fl541) . the luminosity of a 2D black hole is 




'0 


To g {E/rrip) J 27 ih' 


(56) 
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where we use eqn. fl2^ for Tg// and Tq = with n = Note gs is the number of 

polarization which is 1 for scalars and 2 for spin-1/2 fermion and vector bosons. Defining 


E f {E/mp) 
Tog{E/mp)' 


(57) 


we hnd 


Using eqn. 


E = Tnuh 


uTn 


nir 


(58) 


to change variables in eqn. (|56|l . we have for the luminosity 


L = 


9sT^ 

2TTmi 


h 


uTn 


rur. 


h 


uTn 


TUr, 


+ 


uTn 


rUr. 


h' 


uTn 


rUr^ 


un {u) du, 


(59) 


where we define 



(60) 


Note that the change of variables given in eqn. fl57)l is legitimate for the integral if the 
function is monotonic over (0,a;cr), which is the case for the AC dispersion relation. 

Now investigate properties of the luminosity for the AC dispersion relation given in eqn. 
([3]). If To iTip, eqn. fl59l) becomes 

. ... 9sn 


^2 
-0 

27rm^ 


n -|- 2 / uTq 


rrir. 


un (u) du 


(61) 


where we set Umax = oo. For the emission of species of massless particles of spin s, we 
have 


}. (62) 

where Lig {z) is the polylogarithm of order s and argument 2 ; and ({z) is the zeta function. 
li 7] > 0{ri < 0), the luminosity becomes smaller (larger) than that in the usual case with 
g = 0. This is expected from eqn. fl2T)) where Tg// is lowered (raised) due to rainbow gravity 
ii g > 0{g < 0). 

During the late stage of the black hole evaporation process when -C min{l,Xcr}, 

one has = m-Mct ^ hence 

Trip Trip 


3 (ti + 1)! (n + 2) ! Tq 


TT^ 


rUr 


Lin+2 (1) (^0 + 2ui) -|- 


2 - 2 "’ 
u -|- 2 


-C (n 2) ni/2 




nip M — Mcr 
To 


1 + 


g f M - Me, 


nir 


(63) 
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The luminosity of radiation of one species of bosons is 


L = 


gsTo M - Me, 


27imr 


rur. 


1 - 


7 ] 1 

2 n + 1 


M -Mr, 


rUr. 


+ 0 


M -Mr, 


rUr. 


2n' 


(64) 


and that of radiation of one species of fermions is 


dvr 


Qs ( M - Mcr 


rrir. 


1 + 0 


M -Mr. 


TUr, 


2n' 


(65) 


where Qs is the number of polarization. From eqns. fl6T|) and ([65]), we find that the black 
hole evaporates mostly via bosons in the late stage of the black hole evaporation process. 


B. 4D Spherically Symmetric Black Hole 


For a 4D spherically symmetric black hole with hat {x) dx'^dx^ = dO"^ + sin^ 9d(j)‘^ in eqn. 
( 0 . we have found A in eqn. (l2T]l is given by |40| 


A = 



( 66 ) 


where Z = 0,1, • • • is the angular momentum. Since Pr in eqn. fl2T|) is always a real number 
in the geometric optics approximation, one has an upper bound on A 


^ ^ C{r^) 

~ B{r) g‘^{E/mp) 


(67) 


Suppose has a minimum at Team- If the particles overcome the angular momentum 

barrier and get absorbed by the black hole, one has 


^ C {rjj P {Ejmp) 
~ 5(?^min) g^iE/mp) 


( 68 ) 


The luminosity is 

L = (2Z + 1) / 

7 


En 


dE 


= 9s 


EdE 

2TTh^ 


n 


9s C'(rg,iJ 
277^,3 B (rrr,in) 


^'^2nh 
E f{E/mp) 

Tog{E/mp) 

Ef{E/mp] 



( 

/ d 


Jo 

[V 2; J 


n 


P {E/mp) 
To 9 {E/nip) ) g^{E/nip) 


EME, 


(69) 
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where we use eqn. fl5T]) for ue^i- Making change of variables for eqn. (I69|) 


u = 


E f {E/rrip) 


Tog{E/mp)' 

the luminosity of a 4D spherically symmetric black hole becomes 


(70) 


L = 


gsT^ g(^min) r 
27rm^ B (rmin) rnl Jq 


h 


uTn 


rur. 


h 


uTn 


rrir. 


+ 


uTq\ (uT^ 


rUr. 


h' 


nir, 


u^n {u) du. (71) 


In what follows, we focus on the AC dispersion relation. If Tq nip, eqn. GD becomes 


27rm4 B (r^in) rnl Jq 


n + 2 ( uTc 


-T] 


nir. 


vJn iu) du, 


(72) 


where we set Umax = C) 0 . For the emission of species of massless particles of spin s, we 
have 


L^t(2±\ + + 


30 \mp) m25(rmin) 
15 (n + 3)! (n + 2) / Tq 


27r4 


(^) ^ + (2 - 2 C (n + 4) ni/ 2 ] I , 


(73) 


where Lig (z) is the polylogarithm of order s and argument z and ({z) is the zeta function. 
In the late stage of the black hole evaporation process with ^ Xcr}, one hnds 

that the luminosity of radiation of one species of bosons is 


L = 


9sTo C{ri^) fM-M, 


OTrrUp mlB (rmin) V 


1 + 


3g 


2 (n + 3) \ m. 


M -Mr. 


+ 0 


'M - Mr. 





and that of radiation of one species of fermions 

U + 4V mp J y\ nip ) yj 

(75) 

where gs is the number of polarization. Similar to a 2D black hole, the evaporation via 
bosons dominates the 4D spherically symmetric black hole evaporation process in the late 
stage. 

In the geometric optics approximation, a 4D spherically symmetric black hole can be 
described as a black sphere for absorbing particles. The total luminosity are determined 


gg C(rTj 

871 m^B \ nip ) 
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by the radius of the black sphere R and the temperature of the black hole T. For the AC 
dispersion relation, we have for massless particles 


R 



C {rjj 

B (r^in) 1 


1 

r]{E/m^Y 


and T 


Teis = T„^l- ri{E/rn^)". 


(76) 


If T] > 0(?7 < 0), the radius R becomes larger (smaller) than that in the usual case while 
the effective temperature becomes lower(higher) due to rainbow gravity. For the sublu¬ 
minal case with rj > 0(the superluminal case with rj < 0), the competition between the 
increased(decreased) radius and the decreased(increased) temperature determines whether 
the luminosity would increase or decrease. For Tq -C rup, it appears from eqn. d73|) that the 
effects of decreased(increased) temperature wins the competition and hence the luminosity 
tends to become smaller (larger). However for the late stage of the evaporation process, it 
seems from eqns. fl71|l and fITSD that the effects of increased(decreased) radius wins the 
competition and hence the luminosity becomes larger (smaller). 

We now work with a Schwarzschild black hole to investigate more properties of the black 
hole’s luminosity. For a Schwarzschild black hole, one has B (r) = 1 — Vh = 2M, 
K = ^, To = r^in = 3M and = 27M‘^. For M > rup, it shows form 

eqn. (1731) that the corrections to the luminosity from rainbow gravity effects are around 
O Nevertheless, these corrections begin to become appreciable around M < M* = 

\r]\ ' TUp when the second term in the bracket of eqn. fl73|) becomes comparable to 
1. Here c„ is the numerical factor in front of Ug of the second term in the bracket and 
~ 5 - 10 for 1 < n < 10. If = 0, in the late stage of the evaporation process 
with ^ -C min{l,Xcr}, the corrections to the luminosity are around O Similarly, 

these corrections are important when M > M** = nip. Thus, we can conclude that 

for cases with = 0, the effects of rainbow gravity impacts the black hole’s luminosity 
noticeably when M** ^ M < M*. For cases with non-vanishing , the the black hole’s 

luminosity deviates from that in the usual case appreciably when M < M*. In FIG. [6l we 
plot the luminosity L of radiation of one species of bosons against M/nip for examples with 
r; = 0, {r],n) = (1,2), (77,n) = (—1,1), {jq^n) = (—1,2), and {rf^n) = (—1,4). Note that 
the effects of rainbow gravity does not change the black hole’s luminosity appreciably in 
FIG. Elsince M* ~ nip = nip M** for \ri\ = 1. Due to the requirement that 

the energy E of emitted particles could not exceed M — , the luminosities in all cases 

approach zero as M —)■ . 
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FIG. 6: Plot of the luminosity L of radiation of one species of bosons against the mass M/rup 
for a Schwarzschild black hole. The blue line is the usual case. All the other lines asymptotically 
approach the blue line as M S> 1. For the cases with the minimal mass = 0 , the black lines 
asymptotically approach the blue line as M ^ 0. For all the lines, the luminosities approach zero 
as M —)■ 


To study the sensitivity of the BH dynamics to the parameter r], we plot the luminosity 
L of radiation of one species of bosons against M/rUp for different values of rj in FIG. [3 
Since M* ~ \ri\^^^ rrip and M** ~ rup, we hnd that the larger |? 7 | is, the more apparent 

the effects of rainbow gravity on the luminosities become. In FIG. [3 we plot three cases as 
follows: 


(a) Subluminal cases where rj > 0 and = 0. An example with n = 2 is plotted in 


FIG. 7(a) for rj = 1,10,100, and 1000. The red solid line (?7 = 1000) can be barely 


seen since it is too close to the horizontal axis. It is evident that the luminosity for 
M** ^ M < M* gets more suppressed as the parameter r] becomes larger. In other 
words, the effects of rainbow gravity could dramatically slow down the evaporation 
process of the black hole for large enough r]. Hence, the corresponding characteristic 
time scale for the black hole to evaporate from M* to M** could be much longer than 
in the usual case. 

(6) Superluminal cases with the non-vanishing remnant mass Mj.*" where rj < 0 and 


n > 2. An example with u = 4 is plotted in FIG. 7(b) for r] = —1, —10, —100, and 
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(a) Subluminal cases 
(ry > 0) 


(b) Superluminal cases with remnant 
(77 < 0 and n > 2 ) 



(c) Superluminal cases without remnant 
(77 < 0 and 0 < 71 < 2) 


(d) Superluminal cases without remnant 
(77 < 0 and 0 < n < 2) 


FIG. 7: Plots of the luminosity L of radiation of one species of bosons against the mass M/rUp 
for a Schwarzschild black hole. The luminosity L is plotted in subluminal cases and superluminal 
cases with and without remnant for various values of 77 . 


— 1000. In this case, eqn. (1571) gives that ^57*. It shows that the 

luminosity starts to deviate from that in the usual case when M is close to and 
then decreases to zero once M = . Note that becomes larger as we increase 

Ihl- 


(c) Superluminal cases with = 0 where 77 < 0 and 0 < n < 2 . An example with 


n = 1 is plotted in FIGs. 7(c) and 7(d) for rj = — 1 , —10, —100, and —1000. Opposite 
to the subluminal cases, the effects of rainbow gravity could dramatically boost the 
luminosity for M** ^ M < M* in this case if It^I is large enough. For example, the 
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maximum value of the luminosity for r] = —1000 is approximately 10^ times greater 
than that in the usual case and 10^° times greater for r] = —10^®, for which the energy 
scale of Lorentz-invariance violation is assumed around ITeV. The rainbow gravity 
speeds up the process of evaporation and hence the characteristic time scale for the 
black hole to evaporate from M* to M** becomes shorter than in the usual case. 


The luminosity for a Schwarzschild black hole has also been calculated in the framework 
of rainbow gravity in 32| , where massless bosons were considered and the authors used the 
MDR of the form of 

- p‘^ {1 - Xp) = (77) 

’^contrast, there are a number of differences between the calculations in our paper and in 


32| . which are as follows: 


1. The geometric optics approximation have been used to calculate the luminosity of a 
black hole in both papers. In such approximation, a Schwarzschild black hole can be 
described as a black sphere of the radius R and the temperature T. As a result, when 
the luminosity is calculated in the framework of rainbow gravity, we have considered 
effects of rainbow gravity on both R and T. On the other hand, only effects of rainbow 


gravity on T were considered in 


32| . As discussed before, the corrections to R could 


play an important role in the late stage of the evaporation process. This can be 


illustrated by FIG. 7(d), where the corrections to R dominate over those to T and 


make the luminosity smaller than in the usual case for small enough M/mp, although 
the corrections to T tend to increase the luminosity. 

2. In our paper, we used detailed balance condition to show that the average number 
riE^i in the mode with E and i is 

E 


nE,i = n 


T, 


eff 


(78) 


where Te// is given in eqn. (|2T|l . In contrast, the authors of 32|] assumed an aver¬ 


age behavior for particles described 
obtained in section IIIII Therefore in 




a unique average temperature, which is T 


sc 


32l |. the average number ue^z was given by 


nE,i - n I . 


(79) 
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3. We have required t he p articles’ energy E < min {m^Xcr, M — Mcr} while only E < 


mpXcr was used in 


32| . The extra requirement E < M — could dramatically 


change behaviors of the evaporation process in the late stage. In fact, for the usual 
case with rj = 0, Xcr = oo, and M^r = 0, one hnds the luminosity of bosons without 
imposing E < is L ~ which implies that the black hole would evaporate 

completely in hnite time and have a hnal explosion at M = 0. However, if the 
requirement E < M — is imposed, the luminosity L ~ for small enough M. 

Thus, it takes inhnite time for the black hole to evaporate completely and death of 
the black hole is much milder. In our paper, since a black hole in rainbow gravity 
shares the similar the late-stage behaviors with that in the usual case, the lifetime of 
the rainbow black hole is also inhnite. In contrast, the lifetime of a rainbow black hole 
in 32j turned out to be hnite. 


C. 4D Cylindrically Symmetric Black Hole 


For a 4D cylindrically symmetric black hole with hab ( 2 :) dx°‘dx^ = dO^ + a^dz^ in eqn. 
), we have found that A in eqn. fl^ is [40 1 


A = fh^ + 4. 


( 80 ) 


where j is the angular momentum along ; 2 -axis and Jz is a constant. To count the number 
of modes of radiation, we assume the length of the black string is a. Thus, the periodicity 
condition along z-axis gives 

OTrhh 

( 81 ) 


2TTkh , , , ^ 

Jz = -with k & Z. 


In the geometric optics approximation, eqn. fl2lll puts an upper bound on A 


A ^ Ama.-jc — 


_ CjrEj p {E/mp) 

g^{E/mp) 


E^ 


(82) 
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C(r^) 

where has a minimum at rmin- The luminosity per unit length I is 


9s 


l = - = 
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(^9s 
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En 


E f [E/nip) \ dE 


TQg{E/mp) J 2Trh 


ATr'^h^ J 

(^9s C'(rTn) 


d (jh) d ( — 


En 


X « 

^^max -f2 


dE 


f (E/mp) 3 


E f{E/mp) 
Tog{E/mp) 

E f {E/mp 


2TTh^ B (rmin) Jo 5'^ {E/mp) 

^9sn c{ti^) r-\f^ 

27rm^ B (rmin) ml Jq V 


E n 


To g {E/mp) 


dE 


h 


uTn 


mr. 


+ 


mTo 




h' 


uTn 


mr. 


u^n {u) du, (83) 


where u = EiMHEl 
Wneie U Tog{E/mp)- 

For a static uncharged black string 

3am, 


To = 




, one has 


= oo, and 


C {riJ 


= a 


-2 


4vr “ ’ .. 

where b = AM. Since the length of the black string is infinite, one only has E < E^ax = 
mpXcr- The luminosity per unit length for a static uncharged black string is 


I = a 
h 


512%^ 
/ 3am- 


-bs 


0-1 9 9 47Tycr / 

3lm^a^gs^4 f 3cmpb^/^ , f 3amp 




Att 


51/3 


u 


\ Att 




47r 
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3am^ 

Air 


51/3 


u 


u^n {u) du. 


(84) 


In the cases with = 0, if Tq mp {M 1), for the emission of species of massless 
particles of spin s, we have 


I 


81m^a^. 4 


a 


512^5 m + 2ni + ni/ 2 ) 
3 (n + 1)! (n + 2) / 3amp^^/^^ 


TT^ 


\ 47r 


u g 


2-2-^ 

Lin +2 (1) (^^0 + 2ni) H--C {n + 2) ni /2 

n + 2 


( 86 ) 


where Lis {z) is the polylogarithm of order s and argument z and ({z) is the zeta function. 
From eqn. flH^ . we hnd that ^ = I for M 1.Consequently, just like the usual case, 

we have for the rainbow black string that M ~ t~^, which means that its lifetime is inhnite. 
In the following, we will compute the asymptotic value of I when Tq 3> mp {M S> 1). In 

4 

the usual case, the luminosity per unit length I oc Ms. For other cases, the results show as 
follows: 


25 


































Subluminal cases where rj > 0 and = 0. Using 


X = yh{y) r^T] 


-l/n _ V 


— 3ln 


ny^ 


( 86 ) 


in eqn. 


, we have for the luminosity per unit length 


I 


V "9s 
Tmrriia 


n[u, 


u 


du. 


(87) 


^ " 3amp(4M)l/3 

Thus, the luminosity per unit length of radiation of one species of bosons is 


I 


3 r/ "Qs 


(4M) 


1/3 


( 88 ) 


and that of radiation of one species of fermions is 


I 


V *9s 

driTTnipa 


In M. 


(89) 


The luminosity per unit length I is lower than that in the usual case for large M. In 
addition, the radiation of bosons dominates the evaporation process for M 3> 1. It is 
evident that I becomes smaller as r] is increased. An example of one species of bosons 


with rj = 1 and n = 2 is plotted as a black line in FIG. 8(a) Additionally, we plot I 


against M for the examples with n = 4 in FIG. 8(b) for r/ = 1,10,100, and 1000. 


{b) 


Superluminal cases with the non-vanishing remnant mass where rj <0 and n>2 

In this case, one has for radiation of one species of bosons 

3c/, (4M)^/^ 


^ Stt^ rrip Jo 
and radiation of one species of fermions 

9s 


h {y) [h (y) + yh' (?/)] y^dy, 


r*ycr 


I 


271 amj Jo 


h iy) [h {y) + yh' (i/)] y^dy. 


(90) 


(91) 


Similar to the subluminal cases, I could become much lower than that in the usual case 
for large enough M and the radiation of bosons dominates the evaporation process for 
M S> 1. An example of one species of bosons with rj = —1 and n = 4 is plotted as a 


red solid line in FIG. 8(a 


(c) 


Superluminal cases with = 0 where rj < 0 and 0 < n < 2. Using 

1 2 

X = yh (y) ~ (— 77 ) 2^ 7/2^ 


( 92 ) 
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(a) Subluminal cases and superluminal cases with (b) Subluminal cases for various values of r] 
and without remnant 


FIG. 8: Plots of the luminosity per unit length I of radiation of one species of bosons against the 
mass per unit length M for a static uncharged black string. 


in eqn. 


, we find 

^ {-r])^ g^Slmla^ / SarUp 
2 — n 2567 r® 1 47 r 


2 + n 


(4M)3(2->i) / u^n{u)u‘^-^du. (93) 


Since ^ 


^(^ 2 -n) '' 3 for 0 < n < 2 , / could become much larger than that in the usual 
case for large enough M. The larger \r]\ is, the faster the black string evaporates for 
M i:g> 1. An example of one species of bosons with 77 = —1 and n = 1 is plotted as a 


black dashed line in FIG. 8 (a' 


V. CONCLUSION 

In this paper, we have analyzed the effects of rainbow gravity on the temperatures, en¬ 
tropies and luminosities of black holes. Using the Hamilton-Jacobi method for scalars, spin 
1/2 fermions and vector bosons, we hrst obtained the effective temperature T^ff of a black 
hole, which depends on the energy E of emitted particles. By relating the momentum p 
of particles to the event horizon radius Vh of the black hole, the temperatures of a rainbow 
Schwarzschild black hole and a rainbow static uncharged black string were calculated. Fo¬ 
cusing on the AC dispersion relation with f [x) = 1 and g {x) = we computed 

their minimum masses Mcr and hnal temperatures Ter for different values of rj and n. All the 
results were listed in TABLE [U In addition, a non-vanishing minimum mass indicates the 


27 


















existence of the black hole’s remnant, which could shed light on the “information paradox”. 
The entropies were also studied in section IIIII 

In section IIVI we used the geometric optics approximation to compute luminosities of a 
2D black hole, a Schwarzschild one and a static uncharged black string in the framework 
of rainbow gravity. It was found that the luminosity of the rainbow Schwarzschild black 
hole with the mass M deviates from that in the usual case only when M** < M < M* 
for = 0 or < M < M* for a non-vanishing where M* ~ rUp and 

M** ~ nip. In the subluminal cases where p > 0 and = 0, FIG. 7(a) shows that 

the effects of rainbow gravity could signihcantly suppress the luminosity for large r]. In the 


superluminal cases with = 0 where rj < 0 and 0 < n < 2, FIG. 7(c) shows that the 


luminosity could be signihcantly boosted for large \ri\. Similar results for the rainbow static 
uncharged black string with the mass per unit length M S> 1 were also obtained for the 
sub luminal cases and superluminal cases with = 0. 

If the energy scale of Lorentz-invariance violation is A, the naturalness in effective held 
theories implies that \r]\ ~ ^ Schwarzschild black hole, the ehects of rainbow 

gravity start to play an important role when the mass of the black hole M < rup ~ A. 
If Lorentz invariance is violated by quantum gravity, the natural scale A ~ nip. Gurrently, 
there are no experimental evidence that Lorentz symmetry is violated in nature, which 
might suggest that lTeV< A < nip. The experirnental and observational constraints on 


Lorentz-invariance violation are reviewed in 
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5(ll |. As a result, the rainbow gravity plays 


a negligible role for stellar or galactic supermassive black holes. However, the possible 
production of TeV-scale black holes at the LHG or ultra-high-energy cosmic ray collisions 
is predicted by low-scale quantum gravity Mi- Another possible source of small black holes 
is primordial black holes, which are created by primordial density fluctuations in the early 
universe and evaporate for enough long time. Future studies investigating the implications 
of our results on the rich phenomenology of these small black holes would be interesting. 
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